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We consider the single hole dynamics in a modified t — J 
model at finite temperature. The modified model includes 
a next nearest (t') and next-next nearest (t") hopping. The 
model has been considered before in the zero temperature 
limit to explain angle resolved photo-emission measurements. 
We extend this consideration to the case of finite temperature 
where long-range anti-ferromagnetic order is destroyed, us- 
ing the self-consistent Born approximation. The Dyson equa- 
tion which relates the single hole Green's functions for a fixed 
pseudo-spin and for fixed spin is derived. The Green's func- 
tion with fixed pseudo-spin is infrared stable but the Green's 
function with fixed spin is close to an infrared divergency. 
We demonstrate how to renormalize this Green's function in 
order to assure numerical convergence. At non-zero tempera- 
ture the quasi-particle peaks are found to shift down in energy 
and to be broadened. 

PACS numbers: 75.50.Ee, 75.10.Jm, 



I. INTRODUCTION 

Recent angle resolvecLphoto-emission (ARPES) mea- 
surements by Wells et aB and by Pothuizen et aa on the 
insulating Copper Oxide Si^CuO^C^ provide an unique 
experimental probe of the properties of a single hole in an 
anti-ferromagnetic background. Theoretically this prob- 
lem was analyzed in terms of a t — t' — J model .using ex- 
act diagonalization techniques for small clusterso and the 
self-consistent Born approximation (SCBAp. From re- 
cent evaluation of the hopping integrals it was concluded 
that the next-next nearest neighbor hopping matrix ele- 
ment t" is significant and almost as large as the diagonal 
hopping matrix element t', so thatri" should also be in- 
corporated ka model Hamiltonianfl. This was done in a 
recent papem, where also the leading corrections to the 
SCBA were evaluated. nn 

The ARPES experimentsEra are carried out at a tem- 
perature of 300-350K, which is above the Neel temper- 
ature of this compound. Theoretical treatments up to 
now, however, restrict themselves to the zero tempera- 
ture limit, assuming long range anti-ferromagnetic order, 
and spectra are artificially broadened in order to com- 
pare with experiment. It is therefore important to ex- 
tend the SCBA calculations to finite temperature, where 
long range anti-ferromagnetic order is lacking, although 
at room temperature the magnetic correlation length for 
this compound is still about 60 lattice spacings and no 
drastic deviation of the ARPES spectrum at room tem- 
perature from the one at zero temperature is expected. 



Another motivation for this work is that a SCBA tech- 
nique that can cope with the absence of long range mag- 
netic order may be extended in the future to describe the 
spin liquid state of the doped copper oxides. 

We consider a two-dimensional t — t' — t" — J model 
at finite temperature. We apply the modified spin-wave 
theory suggested by Takahashi for 2D Heisenberg model 
at nonzero temperatureQ to deal with a state without 
long range anti-ferromagnetic order. The Hamiltonian 
for t — t' — t" — J model is 



(ij)cr 
J ^ SjSj, 



(1) 



where c\ a is the creation operator of an electron with spin 
a {a =|, J,) at site i on the two-dimensional square lattice, 
(ij) represents nearest neighbor sites, (iji) next nearest 
neighbor sites (diagonal), and (ija) represents next-next 
nearest sites. The spin operator is Si = \c ia a a fjCi^. The 
exchange derived from two magnon Raman scattering is 
J = 1257neyEra. Following the most recent calculation of 
the hopping matrix elements performed by Andersen et 
aE we take: t = 386meV, t' = -105meF, t" = 86meV. 
We set J = 1, so that in these units 



t = 3.1, t' = -0.8, t" = 0.7 



(2) 



We first calculate the hole Green's function with fixed 
pseudo-spin at finite temperature, introduced by a con- 
straint on the sub-lattice magnetization, and evaluate the 
contribution to the self energy due to the virtual absorp- 
tion of spin-waves by the hole. Then the hole Green's 
function with fixed spin, that corresponds to Green's 
function that is actually measured in ARPES, is calcu- 
lated. This Green's function turns out to be close an 
infrared divergency and we show that this instability can 
be avoided by a proper renormalization, assuring that the 
results numerically converge even when a rather limited 
number of grid-points is used. We find that at non-zero 
temperature the quasi-particle peaks broaden and shift 
to lower energy. The shift is independent of momentum 
and is due to the larger effective hole bandwidth as at 
finite temperature the spin order is frustrated. 



II. HOLE GREEN'S FUNCTION WITH FIXED 
PSEUDO-SPIN G D 

At half filling (one electron per site) the model under 
consideration is equivalent to a Heisenberg model. We 
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are interested in the situation when one electron is re- 
moved from the system, when a single hole is produced. 
The dynamics of a single hole in an antLferromagnetic 
background can be described by SCBABEl This ap- 
proximation works very well due to the absencepLa single 
loop corrections to the hole-spin-wave vertexll3ilj. Now 
we have to modify SCBA for finite temperature. The 
main complication is that at finite temperature there is 
no long range anti-Jjerromagnetic order. Nevertheless, 
following Takahashu we introduce artificially two sub- 
lattices: up and down. The bare hole operator di is de- 
fined so that d\ oc Cj-f on the f sub-lattice and oc c$j on 
the I sub-lattice. In momentum representation 



JV(l/2 + m) 



'ST 



<k[ 



N(l/2 + m] 



E c ^ e 



(3) 



where N is number of sites and m — (Si Z ) — is the av- 
erage magnetization. The brackets () denote both quan- 
tum and statistical averaging. The quasi-momentum k 
is restricted to be inside the magnetic Brillouin zone: 
7k = ^ (cos k x + cos k y ) > 0. In this notations it looks like 
dko- has spin a = ±1/2, but actually rotation invariance 
is violated and a is a pseudo-spin which just labels the 
two sub-lattices. Nevertheless the pseudo-spin gives the 
correct value of the spin z-projection: S z = a = ±1/2. 
The coefficients in (||) provide the correct normalization: 



^i4|> = |E( c lT^T>=2^ + <5. 



= 1. 



(4) 



The retarded hole Green's function is defined as 



G d (e,k) = -i (dwa(T)dl a (0))e l ^dT (5) 
Jo 

The t', t" terms in the Hamiltonian (|l|) correspond to the 
hole hopping inside one sub-lattice. This gives the bare 
hole dispersion 



£ok = 4i' cos k x cos k y + 2f"(cos 2k x + cos 2k y 



(6) 



where 7^ = ^(cosk x — cosk y ), /3 i = 4(2t" + t'), and 
(3q2 = 4(2t" — t'). In equation (||) we took into account 
that the sign of a hole dispersion is opposite to that for 
an electron (the maximum of the electron band corre- 
sponds to the minimum of the hole band), and omitted 
the constant term. The bare hole Green's function is 



G d(e, k) = 



1 



e - £ok + i0 



(7) 



For spin excitations we userthe modified spin- wave the- 
oryu (see also review paperta). In order to treat the 



Heisenberg term in the Hamiltonian ([[]) within spin- wave 
theory, it is convenient to use the Dyson-Maleev trans- 
formationta for a localized spin S — 1/2, 

Sf — a], = (2S - ajai)ai, 

S[ — S — ajai, for I £ up sublattice; (8) 

S m = b m , — b m (2S — b m b m ), 

Si^ = —S + b^bm, for m G down sublattice, 

and the Fourier representation for ai and b m : 



(9) 



The summation over q, here and everywhere below, is 
restricted inside the magnetic Brillouin zone. There are 
essentially two ways to find an effective Hamiltonian 
quadratic in the operators a and b. The first way is just to 
drop the quartic terms as is done in linear spin- wave the- 
ory (LSWT). The second way is to treat the quartic terms 
at mean-field level a^ab^b — > (aJa)wb + {(Vw)ab+ cor- 
responding to mean field spin-wave theory. As both ap- 
proximations give very close results we choose use LSWT 
because it is simpler in the present context. 



A. Finite Temperature Correction for Gd 

So far we followed the zero temperature derivation for 
the SCBA. We take the approach of calculating the finite 
temperature corrections to the hole Green's function G^ 
with a diagrammatic, perturbative, method. This frame- 
work can be used if the number of spin- waves per site is 
small, i.e. if T / J is not too large. This is a reasonable 
pre-requisition as in the experiments T / J is of the order 
of 1/4. We will show later on that at these temperatures 
the number of spin-waves is actually small, justifying our 
approach. The advantage of this approach over perform- 
ing the calculations at imaginary, Matsubara, frequencies 
is that in our framework the absence of long-range anti- 
ferromagnetic order is specifically built in. Also problems 
with the numerical analytic continuation on the real fre- 
quency axis are avoided. As is shown later on, an almost- 
divergency occurs that demands large lattices in order to 
get stable results. By performing the calculation on the 
real axis from the beginning, the Green's function can be 
renormalized, lifting the almost-divergency. 

At non-zero temperatures the long range anti- 
ferromagnetic order is destroyed. Very many long- wave- 
length spin waves are excited, and one must take into ac- 
count their non-linear interaction. An approximate way 
to do it is to impose an additional condition that the 
sub-lattice magnetization is zerd3. 
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(\ - a\a x + ~ - blb m ) = 0. (10) 



This constraint gives an effective cutoff of unphysical 
states in Dyson-Maleev transformation. The constraint 
( |io| ) is introduced into the Hamiltonian via a Lagrange 



multiplier \v . Now we must diagonalizc 



(11) 



where A = l + ^ 2 /8. This can be done by the Bogoliubov 
transformation 



fl q = UqQfq + Wq/?L q , 



(12) 



and we find the effective spectrum and Bogoliubov pa- 
rameters 



w„ q = 2, A 2 - 7, 



A 1 



(13) 



In this calculation the usual mean field factorization ap- 
proximation 



> w (a q ctj_Cj T )(c| T Cit) = 1/2 ■ (aqSj ) 



has been used. For simplicity we have omitted in 



the calculation (17) the standard Bose statistics factor 
y/l + riq, but they are certainly taken into account in 
the calculation of diagrams. We stress that the vertex 
( |l7j ) has the same form as in the case of zero temper- 
ature, except for the pseudo-gap in the Bogoliubov pa- 
rameters (|lj). We remind of the fact that at zero tem- 
perature 3k,q=o = because of the Goldstone theorem. 
In the present case due to the pseudo-gap <7k,q=o 7^ 0. 
This is a reflection of the fact that the long range 
anti-ferromagnetic order is destroyed. However, as the 
pseudo-gap is small, its presence will not give rise to large 
effects in the spectra. 




Vq = sign(jq 



These equation show that at non-zero temperature the 
spin-wave spectrum has a gap u\J 1 + ~ v. This elu- 
cidates the meaning of the constraint and the Lagrange 
multiplier. Taking into account that in thermal equilib- 
rium 

Uq = ««q) = (P\Pq) = —fcljT) - 1 

we obtain from (|Io| ) the equation for v 

= 1-|fE— (l + 2"q). 
N ^ UJ v q 



(15) 



This equation gives an exponentially small v, and hence 
an exponentially large magnetic correlation length oc 

Hopping to a nearest neighbor in the Hamiltonian (fil) 
gives an interaction of the hole with spin- waves. 



Hh,sw = y^ffk, q (<4 +a .i.rfkTQ! q 

k,q 

+ 4+qT d k|/3q +H.C.). 



(16) 



In diagrams we will denote the vertex <7k,q by a dot, see 
figure [Bp. This vertex is given by 



ffk.q = (a q d kT |if t |4 +qi ) 



4 *V iv^ 7kUq + 7k+qUq ^ 



(17) 



FIG. 1. Self energy at finite temperature. The solid line 
represents the hole Green's function Gd, and the dashed line 
represents the spin-wave. The first diagram corresponds to 
virtual emission of the spin wave, and the second diagram 
corresponds to virtual absorption. 

Similar to the zero temperature case the spin structure 
of the interaction (|l6|) forbids single loop corrections to 
the hole-spin- wave vertex, so that the self energy is of the 
form 

S ( £ ' k ) = X^ 1 + n q).9k- q ,qGd(e - w q ,k - q) 
q 

+ nq0k, q Gk(e + w q ,k + q)]), (18) 

where first term arises from the virtual emission of the 
spin wave and second term arises from the virtual ab- 
sorption, see figure |[ The self-consistent solution of this 
equation together with the standard relation 



G d (e,k) 



e-e ok -£(e,k)+iO' 



(19) 



gives the retarded Green's function. Due to the definition 
of the operators (||) the Green's function ^ is invariant 
under translation with the inverse vector of the magnetic 
sub-lattice Q = (±ir, ±ir) 



G d {e,k + Q) = G d (e,k) 



(20) 



in spite of the absence of the long range anti- 
ferromagnetic order. 
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B. Results for G^ 

The numerical solution of equation is straightfor- 
ward. To avoid poles we replace iO ~ * iT/2 — i 0.05. The 
energy scale consists of 300 points with variable density 
(concentrated near sharp structures of Go)- The number 
of points in the magnetic Brillouin zone is 10 4 which is 
equivalent to a 140 x 140 lattice. Actually the results are 
almost independent of the grid as soon as it is larger than 
20 x 20. In figure || — i Im G<j(w,k) as a function of u> 
for a cut through the Brillouin zone from k = (7r/2, tt/2) 
to k = (-7T, 0) is shown for two different temperatures. 



< 




(7t,o) — 1 

-7.0 -6.0 -5.0 -4.0 -3.0 
CO 

FIG. 2. Plots of — i Im Gd{u>, k) as a function of uj for a cut 
through the Brillouin zone from k = (n/2,n/2) to k = (tv, 0) 
for T=0.01 and T=0.3. 

The number of spin-waves per site for the highest 
temperature, T = 0.4, we considered, is sw 0.05, so 
much smaller than unity, justifying the perturbational 
approach. We recall that we use the set of parameter* 
(||) based on Reffl. The same set has been used in Ref.H 
for a zero temperature calculation. Quasiparticle ener- 
gies and residues obtained here are quite similar to that 
at zero temperatureQ. We see from figure || that tempera- 
ture shifts the quasi-peaks positions to lower energy and 
broadens them. The explanation for this is that there 
are two contributions to the self-energy at non-zero tem- 
peratures. One term originates from virtual spin- wave 
emission processes and the other from virtual spin- wave 
absorption. The former contribution is also present at 
T=0 and is enhanced at finite temperatures by a factor 
1 + n q , i.e. at finite temperature the interaction between 
the spin- waves and the hole is effectively increased due to 
this process. The matrix element for emission of an extra 
spin-wave is larger if the number of spin-waves already 
present is larger, because of the bosonic nature of the 
spin-waves. The part of the self-energy that is due to the 
spin-wave emission is multiplied by a factor larger than 
one at finite temperature, causing a shift of the quasi- 
particle peaks to lower energy that is nearly uniform in 



the Brillouin zone and a shift of the incoherent part of 
the Greens' function to higher energies. This has a simple 
physical reason. At non-zero temperature the hole prop- 
agates more easily because the anti-ferromagnetic order 
is frustrated. This causes a uniform shift of all quasi- 
particle poles to lower energy and does not effect their 
dispersion, as the quasiparticle dispersion is determined 
predominantly by the magnetic interaction. So the effect 
of non-zero temperature is qualitatively different from 
the effects of doping, where it is found that a reconstruc- 
tion of the quasiparticle dispersion takes place, attributed 
to the frustrated magnetic order in the doped systemllZl. 
The contribution to the self-energy at finite temperature 
because of the spin- wave absorption is mainly responsible 
for the broadening of the quasi-particle peaks. 

III. HOLE GREEN'S FUNCTION WITH FIXED 
SPIN G c 

The operators d^, c?ki discussed in the previous sec- 
tion are defined at different sub-lattices. The operators 
on two sub-lattices are useful as mathematical construc- 
tions but when a photon removes an electron from the 
system, it does not differentiate between the sub-lattices, 
and moreover, at nonzero temperature there are no sub- 
lattices at all. Therefor we have to define the particle 
operator that is relevant for photo-emission independent 
of the sub-lattice. This particle operator is: 

Ck - = \/^E c - e4krs - ( 21 ) 

i 

The normalization is chosen in such a way that 

(0|c kTCkT |0) = |(0|^4 TCiT |Q) = l. (22) 

i 

The corresponding retarded Green's function is 

roc 

G c (e,k) = -i (c k >)c kCT (0)) e * eT dr. (23) 
Jo 

This is the Green's function measured in ARPES. 

Now we have to find the relation between G c (e, k) and 
Gd(e,k). The operator Ck<j acting on half filled ground- 
state produces a single hole. We denote the correspond- 
ing amplitude by dk and denote it in figure ^|a as a cross. 
The thick line corresponds to the Green's function G c 
(^3| ) and the thin line corresponds to the Ga @. The 
amplitude dk equals 



ak = (dkTCki) = 




The operator Ck CT acting on a state of the system can 
also produce a hole + spin-wave. This amplitude is 
shown in figure |^b as a circled cross with the dashed 
line being a spin-wave. We denote this amplitude by 6k, q 
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frk-q,q — (/5qdk-qlCkl 
2^2 



N 
2\/2 



E s . +e ' qri 



(25) 



(/3 q ^(1 - aja^e 




There is some ambiguity in the last step of this deriva- 
tion. If we neglect the quartic term /3 q (l— a|ai)ai — ► /3 q ai 
we get a value of 6k, q by a factor y2 larger than that 
given by eq. ( |25| ) . If we treat the quartic term on a mean 
field level then a\a,i — > 1/2 and we get a value of 6k, q , a 
factor \/2 smaller than in eq. (|25|) . The correct value is 
somewhere in between. We choose the vertex 6k. q to be 
the same as in the zero temperature caseo. We will see 
that this provides the correct sum rule for the Green's 
function G c , and this is a justification of our choice. We 
stress that ( p5| ) is a bare vertex. It corresponds to the 
instantaneous creation of a hole + spin- wave, but not the 
creation of a hole with a subsequent decay into a hole + 
spin- wave. 



A. Finite Temperature Correction for G c 



k-q 



k-tf - 



k-q 



FIG. 3. The vertices: a) - single hole creation, b) - hole 
+ spin-wave creation, c) - hole creation with spin wave an- 
nihilation, g) - usual hole-spin- wave vertex. The thick line 
correspond to G c , and the thin solid line corresponds to Gd- 
The dashed line is the spin-wave. 

At finite temperature there is another possibility: the 
creation of a hole with the absorption of a spin wave from 
the thermal bath. We denote this amplitude by Ck. q - It 
is shown in figure |^c, and for simplicity we also denote it 
as a circled cross. The derivation of Ck, q is quite similar 
to (Ea) and the result is 



Ck,q = («!. 



Jk-qjTki; « \/ JjUq 



(26) 



In eqs. ( |25| ) and ( |26|) we h ave omitted the standard Bose 
statistics factors yj\ + n q in (^J) an d y / ^q in ( pof ) . These 
factors are taken into account separately in the calcula- 
tion of the diagrams. 



+ -8- 



FIG. 4. Dyson equation relating Green's functions G c 
(thick solid line) and Gd (thin solid line). 

Now we can find the relation between the Green's func- 
tions G c ( |23| ) and Gd (||). In leading t approximation it 
is given by the diagrams presented in figure ^ where the 
thin solid line represents the bare hole Green's function 
God (0) . Each self-energy insertion should be understood 
as a combination of spin- wave emission and spin- wave ab- 
sorption diagrams, similar to that in figure |l|. Now let 
us dress these diagrams by higher orders in the hopping 
t. As we already discussed above, there is no single loop 
correction to the "dot" . We neglect double loop correc- 
tions to the "dot" as it has been done in SCBA. Therefore 
the only possibility is an introduction of self energy cor- 
rections to Gd- To take into account all these corrections 
we need just to replace all bare hole Green's functions 
God by dressed hole Green's functions Gd- So, [| actually 
represents a Dyson equation relating G c ( p3| ) and Gd (j^)- 
In analytical form it is 

G c (e,k) = a k G d (e,k)+Si(e,k) 

+ 2a k G d (e, k)S 2 (e, k) + G d (e, k)E§(e, k), (27) 

where 

Si(e, k) = £[(1 + n q )6 k _ q>q G d (e - w q , k - q) 



n q c kiq G d (e + w q ,k + q)], 



(28) 



and 



S 2 (e, k) = yj[(l + ™ q )6k-q,q.9k-q, q G d (e - Uq, k - q) 
q 

+ « q Ck :q 3k, q G d (e + cj q ,k + q)]. (29) 



B. Sum rules 

Let us check now the sum rules. All singularities of the 
retarded Green's functions are in the lower half plane of 
complex e. Therefore if we integrate eq.([l9|) over e from 
—00 to +00, this integral can be replaced by the integral 
over an infinite semi-circle in the upper e half plane. For 
infinite e, Gd = God, and we get the well known sum rule 



1 f 00 

-Im G d e,k)de = l, 

7T J-00 



(30) 



which agrees with eq.(|4|). If we integrate eq.fl27j) in the 
same limits, the terms which contain more than one 
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Green's function give no contribution, because the in- 
tegral can be transferred into the upper complex e half 
plane, and we find 

--Im f G c (e,k)de= (--Im [ G d (e,k)dt 

I" J-oo \ f J 

^4 + K 1 + n q) & k-q,q + ™qCk,q]^ 



N 
q 

l^A ( i + 2nq) = i ! 

N ^ LO va 



(31) 



a lattice of at least 200 x 200, and for lower temperatures 
an even a bigger lattice. 

We can avoid this problem by renormalizing Si, so that 
we can work with a reasonable grid-size. Let us rewrite 
eq. (p8[) in the form 



Si(e,k) = S fl (e,k) + T, RR (e,k), 



(33) 



where 



Sfl(e,k) = £(l + n q )6 2 _ c 



[G rf (e-w q ,k-q)-G rf (e ! k)] (34) 
«q4, q [Gd(e + w q , k + q) - G d (e, k)] , 



where we have used eqs. (|_3() and (|15|). Thus equation 
(p7|) reproduces the correct normalization: (Olc^Ckt |0) = 
1. This also proves that the vertices (^5|) and (^6|) are 
correct. 



The vertices &k,q fl25[ ) and Ck, q (|26|) are invariant un- 
der translation by the inverse vector of magnetic sub- 
lattice Q = (±7r, ±7r): 6k+Q,q = f'k.q- At the same 
time the vertex c/k,q (£JJj) changes sign with this transla- 
tion: <?k+Q,q = — .9k, q- Therefore the self energy £2(5, k) 
changes sign at k — > k + Q and 



G c (e,k + Q) ^G c (e,k). 



(32) 



The imaginary part of G c (e, k) gives directly the spectra 
measured in ARPES experiments. This Green's function 
can be calculated using the Dyson equation ( |27| ) as soon 
as we have found G d in SCBA (pi). 



and 



S^(e,k) = G rf (e,k)^[(l + n q )& 



k-q,q + n q C k,q] 



(35) 



Numerical calculation of S# does not cause any trouble 
because it is well convergent at small q. On the other 
hand Sjyj can be easily calculated analytically using 
the modified spin- wave theory equations ([l3|) and (|15|): 
£#ii(e, k) = iGd(e, k). Using this procedure the calcu- 
lation can be done at arbitrary small temperature. The 
results are practically independent of the grid as soon as 
it is larger than 20 x 20. The plots of —- Im G c (e,k) 
as a functions of e for k = (it/2,itJ2), k = (7r/2,0), and 
k = (tt, 0) are presented in figure \ 



C. Results for G c 

Numerical evaluation of G c at finite temperature, how- 
ever, is more complicated than at T=0. The problem lies 
in the infrared divergence of the integrand of Si(e, k) at 
small q. To clarify this we compare the small q behavior 
of the integrands of the self-energies S, Si, and S2. Small 
q means that 1/£m <C q -C T, where £m °c exp(l.l/T) 
is the magnetic correlation length. In this region the 
spin- wave mean occupation number is n q ~ T/q, and 
the vertices are <7k, q ~ ^/q, &k, q ~ Ck, q ~ l/V?- The 
self energy S has an integrand oc n q g£ q rf 2 9 ~ Tqdq. It 
is convergent at small q and therefore numerical calcu- 
lation of S is straightforward and the finite temperature 
generalization of SCBA is as simple as zero temperature 
SCBA. For the integrand of the self-energy S2 on finds 
oc « q 6k. q .9k,q^ 2( ? ~ Tdq, which is also convergent at small 
q. The situation is different in self-energy Si: it is loga- 
rithmically divergent at small q: oc q <?? ~ Tdq/q. 
There is no real divergence, however, because the inte- 
gral is convergent at q ~ 1/£m, but to calculate this in- 
tegral numerically by "brute force" one needs a grid with 
Aq <C 1/Cm- One then needs for example, when T=0.25, 



/\ *10 
(71,0) //T^^. 








(71/2,71/2) 



-8.0 -6.0 -4.0 -2.0 

CO 



FIG. 5. Plots of — ^ Im G c (lu, k) for different temperatures 
and k-points. The thin line is the result for T = 0.01, the 
thick lines for T — 0.2, 0.3 and 0.4. The spectral weight at 
the (n, 0) point is multiplied by a factor 10. 

The half widths of quasiparticle peaks of G c are slightly 
larger than the half widths of quasiparticle peaks of G d . 
The spectra at T=fU)l quite well agrees with zero tem- 
perature calculation^. We stress that the agreement is 



G 



not trivial. At T=0 long range anti-ferromagnetic order 
is assumed, and in the present work we used a quite differ- 
ent approach based on a state without long range order. 
The agreement indicates that these two approaches are 
consistent. For non-zero temperatures the trend in the 
spectra for G c is the same as for Gd', the temperature ef- 
fect is to shift the quasi-particle peaks uniformly to lower 
energies and to broaden them. The spectra presented in 
figure H should -be directly compared with ARPES ex- 
perimental dataau. They reasonably reproduce positions 
and residues of experimental peaks, but fail to reproduce 
widths of the peaks. 

IV. CONCLUSIONS 

We considered the two-dimensional t—t' — t" — J model 
at finite temperature, and developed a technique to deal 
with the state without long range anti-ferromagnetic or- 
der. There is hope to extend this technique to the spin 
liquid state of the doped system. We generalized the 
self-consistent Born approximation to the case of nonzero 
temperature and derived the Dyson equation which re- 
lates the single hole Green's function with fixed spin to 
the single hole Green's function with fixed pseudo-spin. 
This equation is sensitive to very large distances of order 
of magnetic correlation length and therefore not conve- 
nient for computations. To overcome this problem we 
developed a renormalization procedure which allows one 
to exclude large distances and to work with a relatively 
small lattice: the results are independent of the lattice 
size as soon as it is larger than 20 x 20. The effect of 
a finite temperature is a broadening and a shift of the 
quasi-particle peaks to lower energy, independent of the 
momentum. This is attributed to the frustrated mag- 
netic order at finite temperature. The calculated ARPES 
spectra demonstrate that temperature broadening is not 
enough to explain widths of the experimental spectratm. 
This strengthens the argument that other degrees of free- 
dom contribute to the peak widthEI. 
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